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TECHNICAL NOTE 2365 

ANALYTICAL EVALUATION OF AERODYNAMIC CHARACTERISTICS OF TURBINES 

WITH NONTWISTED ROTOR BLADES 
By William R. Slivka and David H. Silvern 


SUMMARY 

An analytical investigation was made of the aerodynamic char- 
acteristics (relative rotor -entrance Mach number, reaction factor, 
rotor turning angle, radial variation of specific work output, and 
rotor-exit tangential-velocity head) of single-stage turbines with 
nontwisted rotor blades, and an evaluation of these characteristics 
on the basis of the simplif ied-radial-equilibrium analysis was made 
by comparison with the corresponding characteristics of free -vortex 
turbines. In order to obtain an indication of the effects of the 
radial shifts of weight flow on the significant parameters discussed 
in the comparison, further analyses were made of a contemporary single- 
stage turbine where the forces due to the curvature of the streamlines 
in the axial -radial plane are approximated. 

On the basis of the simplified-radial-equilibrium analysis and the 
analysis accounting for streamline curvature in the axial-radial plane, 
the aerodynamic characteristics of nontwisted-rotor -blade turbines are 
approximately those of free -vortex turbines intended for similar appli- 
cation for values of hub -tip -radius ratios that are used in current 
turbines . 


INTRODUCTION 

Nontwisted turbine blades were first used in early steam-turbine 
practice without regard to radial variations in flow conditions. 
Although more recently steam- turbine designers have recognized that 
certain radial variations of flow do exist, consideration of manufac- 
turing costs have generally limited the extent to which these radial 
variations have been accounted for. For aircraft gas turbines, the 
great importance of high single-stage power output, together with 
high turbine efficiency, has required that allowances for radial 
variations in flow conditions be made in design techniques. These 
techniques are generally based on either free -vortex flow or constant 
specific weight flow, either of which results in twisted turbine rotor 
blades. 
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Recent research has indicated that it is desirable to incorporate 

multiple cooling passages within turbine rotor blades (reference 1) . 

The manufacture of such cooled blades will be simplified if the rotor 
blades can be of uniform camber and without twist along the blade 
length, that is, nontwisted rotor blades. 

The design of nontwisted rotor blades in combination with non- 
twisted stator blades is discussed in reference 2. Analysis of the 
data of reference 2 shows that the relative flow angle at the rotor 
entrance varies along the radius. Variation of the relative rotor- 
entrance flow angle, combined with nontwisted rotor blades, however, 
will very likely result in increased losses in a turbine that has high 
aerodynamic loads and a high relative rotor -entrance Mach number. 

In the analysis made at the MCA Lewis laboratory and presented 
herein, the use of nontwisted rotor blades in combination with twisted 
stator blades is considered as a method of obtaining structural 
simplicity and efficient aerodynamic performance of turbines. The 
manufacture of twisted and cooled stator blades is presumably easier 
than twisted and cooled rotor blades because of the freedom in select- 
ing materials and manufacturing processes that the low stator-blade 
stress provides. 

A general method of design analysis for various types of radial- 
flow variation is presented in reference 3 for (a) the simplified- 
radial-equilibrium case where the forces due to the curvature of the 
streamlines in a axial-radial plane are neglected, and (b) the radial- 
equilibrium case where these forces are included. No detailed analysis 
is made, however, for the particular case of nontwisted-rotor -blade 
turbines . 

The purpose of the present report is to investigate the aerodynamic 
characteristics of nontwisted-rotor -blade turbines and to evaluate the 
characteristics of this type of design by comparison with the character- 
istics of free-vortex designs. The free-vortex turbine was chosen for 
comparison because it is one of the most common types used in current 
design practice. 

The evaluation is based on the simplified-radial-equilibrium 
analysis, where the parameters of weight flow, relative rotor -entrance 
Mach number, reaction factor, rotor turning angle, specific work output, 
and rotor-exit tangential -velocity head of single-stage nontwisted- 
rotor -blade turbines are compared with the corresponding parameters of 
single-stage free-vortex turbines. Further analysis is then made of a 
particular single-stage turbine where the forces due to the curvature 


2117 


2117 


NACA TW 2365 


3 


of the streamlines in the axial-radial plane are approximated in order 
to obtain an indication of the effects of the radial shifts of weight 
flow on the parameters investigated in the comparison. 


ANALYSIS 

The aerodynamic design of a turbine may be divided into two phases. 
The first is the analysis of the radial variations of the velocity 
vectors at the stator and rotor exit and the second is the design of 
the blade profiles to achieve the desired velocity diagrams. The 
design of the blade profiles and the effect on the conditions of the 
gas at the rotor entrance and exit are not considered in this analysis. 
Throughout the analysis, the relative air angles at the rotor entrance 
and exit are assumed equal to the rotor -blade angles. The symbols 
used herein are defined in appendix A. 


Simplified-Radial-Equilibrium Analysis 


Flow at rotor entrance. - An analysis of the flow conditions at 
the rotor entrance, with the specification that the relative rotor- 
entrance angle is constant along the radius, is presented in appendix B. 
In this analysis, the following assumptions were made: axial symmetry, 

constant total temperature along the radius, constant entropy along the 
radius, and simplified radial equilibrium. The following equation 
relating the radius ratio r/r Q to the velocity- vector -diagram terms 
(fig. l) is obtained; 




sin 2p 

\[2 


arc tan 


tan ft 

'/2 




(i) 


This equation was obtained by integrating equation (B13) between the 
limits of 1 and V-q/u in the right member and 1 and r/r G in the 
left member. 


Equation (l) was used in the construction of figure 2, where V u /U 
is plotted against r/r Q for constant values of ft. Lines of constant 
stator -exit angle a varying from 10° to 45° are also plotted. The a 
lines were determined from the expression 



tan a = 1 


tan ft 


( 2 ) 
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Figure 2 was used in the evaluation of the nontwisted-rotor -blade 
turbines to determine the radial variations of velocity vectors. The 
chart was used as follows: Values of 0 m , ( V u /u ) m , and U m were 

assigned and the corresponding value of (r/r 0 ) m was then determined. 
Other values of r/r Q were determined from the following expression: 

_r /_r_\ _r_ 

r o \ r o/m r m 

For each value of r/r Q , the corresponding values of Vy/u are 

obtained from the chart for the given value of (3. With these param- 
eters known, the velocity -vector diagrams are then calculated and, with 
given values of stagnation temperature and pressure, the gas properties 
along the radius may be easily calculated. 


In order to determine quickly the aerodynamic characteristics of 
nontwisted-rotor -blade turbines when a number of possible turbine- 
design conditions are being investigated, figure 3 was constructed for 
values of stator-exit angle a from 10° to 45°. A weight-flow param- 
eter - ^ - T — — is plotted against radius ratio — for constant 

P « r o 2 
values of p , 

is defined as 


y = 1.3, and values of M Q from 1.0 to 0.5, where Mq 


u 0 A 


•cr * 


Lines of constant relative rotor-entrance Mach 


number My and a line that divides the charts into supersonic and 
subsonic stator-exit conditions are also plotted. The method of con- 
struction of figure 3 is presented in appendix C. 


The use of figure 3 in conjunction with the use of figure 2 may be 
illustrated as follows: For given values of V u , p, and U at some 
radius, the corresponding value of r/r Q i- s obtained from figure 2. 
With U and r/r Q known, U Q is calculated and, with a given design 

value of T', Mq is then determined. With p' known and with an 
assumed value of flow coefficient that accounts for the effects of 
annulus-wall boundary layer and blade wakes in reducing the flow area, 
the limits of r/r Q that pass the desired weight flow are obtained 
from figure 3. At these limits of r/r Q , the relative rotor-entrance 
Mach number and the stator -exit angle are obtained from figure 3. 


The flow conditions for the free -vortex turbines were determined 
in accordance with the methods presented in reference 3. 
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Flow at rotor exit, - An analysis of the flow conditions at the 
rotor exit, with the specification that the relative rotor-exit angle is 
constant along the radius, is presented in appendix D. In this analy- 
sis, the following assumptions were made: axial symmetry, constant 

entropy along the radius, simplified radial equilibrium, and streamlines 
at the rotor entrance and exit lying on the same radius. The following 
equation defining the change of tangential velocity with radius as 
derived in appendix D applies: 


The conditions at the rotor exit were determined from equation (3) 
in the following manner- The value of was assigned at the mean 

radius and the value of necessary to satisfy continuity was deter- 

mined from figure 4 and the method presented in appendix E. This method 
determines P 3 from an evaluation of flow conditions at the mean radius 
by the assumption that the specific weight flow at the mean radius is 
equal to the average specific weight flow. When the values of 
V u 3 m , P 3 , and rotor -entrance conditions are known, a solution of 

equation (3) was obtained by using the Runge-Kutta method of numerical 
integration, as presented in reference 4. With this information, the 
specific weight flow was then integrated to determine the accuracy of 
the value of P 3 obtained from figure 4. In all cases considered, this 
first approximation was sufficiently accurate. It should be noted that, 
because of the small values of exit tangential velocity existing in 
most cases, the use of large increments of radius yielded sufficiently 
accurate results in the numerical integration. 

The flow conditions for the free -vortex turbines were determined 
in accordance with the methods presented in reference 3. 

Radial Equilibrium Accounting for Streamline Curvature 

In reference 3, it is shown that the assumption of simplified 
radial equilibrium leads to errors in the determination of the radial 
variations of flow. This result especially applies to nontwisted- 
rotor-blade turbines where large radial shifts of weight flow occur. 

This section will be devoted to a method of approximating the forces 
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due to the curvature of the streamlines in the axial-radial plane in 
order to obtain an indication of the effects of the radial shifts of 
weight flow on the results obtained by the simplified-radial-equilibrium 
analysis . 


Flow at rotor entrance. - The equations describing the radial var- 
iation of velocity at the rotor entrance midway between the stator and 
rotor blades for nontwisted-rotor -blade turbines and free -vortex tur- 
bines that include the additional force term due to the curvature of the 
streamlines in the axial-radial plane are presented in appendix F . 

The assumptions used in this analysis were: axial symmetry, constant 

total temperature along the radius, and constant entropy along the 
radius. The equations that define the variation of the axial velocity 
with radius for the cases considered are as follows: 


Nontwisted rotor blade: 
V. 


J ’ x 

d Vu 

d' — 
\ r m 


_1_ 

r 


2+4 


ir) cot p Aw 


2 r 




(4) 


* x P 

— CSC P + cot P 


Free -vortex rotor blade: 


c^x 

dr 



(5) 


where C is the radius of curvature of the streamline in the axial- 
radial plane at the entrance station (fig. l) • It was found to he more 
convenient to present these equations in terms of the radial variation 
of axial velocity instead of the radial variation of tangential 
velocity. 


A simplified method for the solution of these equations for a 
multistage machine is presented in reference 3 with the assumption that 
the form of the streamlines in the axial-radial plane is sinusoidual. 
For a single-stage machine ^ it is possible that the streamlines will he 
more closely approximated hy the equation 


= r 

X^oo 



r 


+ Ae 


2117 


2117 


UACA TTJ 2365 


7 


because this equation will satisfy the condition of the streamlines 
lying on the same radius at ±oo and having a maximum deflection between 
the blade rows. 

It is assumed in this analysis that the maximum deflection of the 
streamlines occurs midway between the stator and the rotor inasmuch as 
it is probable that at this station the specific weight flow has its 
greatest variation. It is further assumed that the inflection point of 
the streamlines occurs at the midchord point of the stator (fig. l) . 

When these assumptions are used, the curvature of the streamlines 
at the point of the streamlines midway between the stator and rotor 
blades equals -A/b 2 , where A is the maximum deflection ( r max - r ) 

and b is the axial distance from the point of inflection to the 
point of maximum deflection. The following equations, which were 
solved by the numerical methods previously used, were obtained by sub- 
stituting for the radius of curvature in equations (4) and (5): 


Nontwisted rotor blades: 



( 6 ) 

Free -vortex rotor blades: 


dV^ -V^ 
dr b 2 


(?) 


Inasmuch as A is not directly determinable, a trial solution is nec- 
essary. The first approximation is found by integrating equations (6) 
and (7) with A equal to zero (simplified-radial-equilibrium case). 
From the velocities and the densities so determined, the weight flow 
was plotted as a function of radius. The radii that divide the weight 
flow into equal parts were determined and the deviation of these radii 
from those determined by dividing the annulus area in equal parts 
(constant specific weight flow at rotor entrance) were found. These 
deviations A were then used to calculate a second approximation. 

This procedure was repeated until the assumed values of A checked 
the final calculated values. After continuity was satisfied, the value 
of V„ ™ did not change with successive approximations of A. 

J\. y ill 
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Flow at rotor exit. - At the rotor exit, it was assumed that the 
effect of the curvature of the streamlines in the axial-radial plane 
is negligible. The analysis as presented in appendix D is therefore 
applicable. However, U 2 is not assumed equal to U 3 as in the 
simplified-radial-equilibrium analysis and equation (D6) was therefore 
used, rather than equation (3) . The conditions at the rotor exit were 
determined from equation (D6) by the same method of solution as was 
previously presented in the section on the flow at the rotor exit for 
the simplified-radial-equilibrium analysis. 


COMPARISON OF RESULTS 
S implif ied -Radial -E qui librium Analy s i s 

The aerodynamic characteristics of weight flow, relative rotor- 
entrance Mach number, reaction factor, rotor turning angle, work out- 
put, and specific rotor-exit tangential-velocity head were determined 
for nontwisted-rotor -blade and free-vortex turbines with the following 
operating conditions; these conditions were chosen to cover the range 
of current design limitations: 


Case 

Stagnation 

inlet 

tempera- 

ture, 

T i 

(°R) 

Mean -radius 
stagnation- 
pressure 
ratio, 

Mean- 

radius 

blade 

speed, 

U m 

(ft/ sec) 

Mean- 
radius 
rotor -exit 
tangential 
velocity, 

(v ,) 

v u,3 m 
(ft/sec) 

Mean-radius 
relative 
rotor - 
entrance 
angle, 

Pm 

m 

(deg) 

Turbine 

internal 

efficiency, 

i 

2500 

1.5 

900 

0 

40, 50, 60 

0.85 

2 

2500 

2.0 

1100 

0 

40, 50, 60 

.85 

3 

2500 

2.5 

1300 

0 

40, 50, 60 

.65 

4 

2500 

3.0 

1400 

0 

40, 50, 60 

.85 

5 

2500 

2.5 

1300 

-0.2 U m 

40, 50, 60 

.85 


A weight-flow ratio w defined as 


P* 

2n L pV r dr 
y/ r h x 
/ \ / 2 2 , 

(P V x>m ” (r t ' r h > 
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is plotted in figure 5 against hub -tip -radius ratio for the preceding 
cases. The significance of the weight-flow ratio is two-fold: (l) It 

demonstrates the deviation of the actual weight flow as compared with 
the weight flow determined from the mean-radius condition; and (2) it 
serves to demonstrate the comparison of the actual weight flow of a 
nontwisted-rotor -blade and a free -vortex turbine, with both turbines 
having the same mean -radius conditions. For all cases considered, 
w = 1 for the free-vortex turbines. For the nontwisted-rotor -blade 
turbines, however, w < 1. This weight-flow ratio decreased as the 
hub -tip -radius ratio decreased, giving a minimum value for the cases 
considered of approximately 0.89 at a hub -tip -radius ratio of 0.6 for 
case 5. At a hub -tip -radius ratio of 0.75, a value that more closely 
approximates the values used in current turbines, the minimum value 
of the weight-flow ratio was approximately 0.96. 

In order to make the comparison of the aerodynamic character- 
istics of the nontwisted-rotor-blade and free-vortex turbines valid 
for equal weight flow for the same hub- tip- radius ratio and mean 
radius, the axial velocity of the free-vortex turbines was decreased so 
that the weight flow of the two types of turbine were equal for 
every hub -tip -radius ratio. This equality resulted in a change in 
the relative rotor -entrance angle at the mean radius for the free- 
vortex turbine at every hub -tip -radius ratio and the remainder of the 
comparison is made on this basis. These results are presented in fig- 
ures 6 to 11, where the parameters of relative rotor -entrance Mach 
number, reaction factor, radial variation of work output, rotor 
turning angle, and specific rotor-exit tangential -velocity head are 
plotted for both the nontwisted-rotor-blade and free-vortex turbines. 
Inasmuch as the magnitude of the values of these parameters for 
cases 2 and 3 were within the limits of the values of the parameters 
for cases 1 and 4, cases 2 and 3 are omitted for the remainder of 
the comparison. 

Relative rotor-entrance Mach number. - The value of the maximum 
allowable relative rotor -entrance Mach number is one of the limiting 
factors in the design of turbines. Because the maximum value of this 
Mach number always occurs at the hub of the machine, the comparison 
of Mach number for the nontwisted-rotor-blade and the free-vortex 
turbines was restricted to comparison at the hub. The relative rotor- 
entrance Mach number at the hub is plotted in figure 6 against hub- 
“tip-radius ratio for the cases considered. For all the cases con- 
sidered, the Mach number for the nontwisted-rotor-blade turbines was 
greater than the Mach number for the free-vortex turbines. At a hub- 
tip-radius ratio of 0.6, the Mach number for the nontwisted-rotor- 
blade turbines was approximately 0.02 to 0.16 greater than the Mach 
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number for the free-vortex turbines; the higher values occur at the 
higher turbine pressure ratios and higher relative rotor- entrance 
angles. The magnitude of the effect on turbine losses associated 
with increases of Mach number of the order discussed is unknown and 
will depend on the blade loading and the critical inlet Mach number 
of the particular blade sections. 

Reaction factor. - The reaction factor t defined by 
= (W 32 -W 22 )/W 2 ^ is a measure of the static-pressure change across 
the rotor. The value of i r based on isentropic flow is equal to 
zero for impulse blading and this value is often considered a limiting 
point in contemporary turbine design. In general, the values of t 
decrease from tip to hub, which makes the hub of the machine critical 
with regard to \|r. The comparison of the reaction factor was there- 
fore restricted to the hub of the machine. A plot of hub reaction 
factor \(r h against hub-tip-radius ratio is shown in figure 7 for the 

nontwisted-rotor-blade and the free-vortex turbines. For cases 1 
and 4, the reaction factors for both types of turbine are approximately 
equal for all hub-tip-radius ratios and relative rotor- entrance angles. 
For case 5, the reaction factor for the free-vortex turbines is from 0 
to approximately 0.2 greater than the reaction factor for the nontwisted- 
rotor-blade turbines. The magnitude of these differences is considered 
negligible. 

Work output. - Two problems are associated with radial variations 
of specific work output: (1) the problem of possible losses associated 

with the mixing to constant conditions and the possible induced losses 
arising from secondary flows associated with radial gradients of work 
output, and (2) the problem of obtaining a given average specific work 
output. The specific work output for the nontwisted-rotor-blade tur- 
bines varies along the radius, whereas the specific work output is 
constant along the radius for the free-vortex turbines, as shown in 
figure 8, where a ratio of specific work to specific work at the mean 
radius is plotted against radius ratio r/r m for the cases considered. 

Cases 1, 4, and 5 show similar variations. For hub-tip-radius ratios 
of 0.6, this variation from the value at the mean radius was of the 
order of 5, 10, and 15 percent for relative rotor-entrance angles of 

40°, 50°, and 60°, respectively. For hub-tip-radius ratios of 0.75, 
which more closely approximate the values used in current turbines, 
this variation is of the order of 2.5, 5, and 7.5 percent for relative 
rotor-entrance angles of 40°, 50°, and 60°, respectively. 


2117 


2117 


NACA 'M 2365 


11 


It is estimated that the mixing losses associated with the varia- 
tion of specific work output presented are negligible, which is based on 
the assumption that only that part of the velocity head that remains 
after the flow is mixed to constant conditions is considered recoverable. 
The induced secondary losses are also estimated to be small. 

In order to determine the effect of the variation of specific work 
output on the over-all work output, an average-specific-work ratio h 
defined as 


h = 



u < v u,2- v u,3>e v x r 
[ u < v u,2- v u,3i] pV ar 

in 


is plotted in figure 9 against hub -tip -radius ratio for the nontwisted- 
rotor -blade and the free-vortex turbines. It may be seen from figure 9 
that, if the nontwist ed-rotor -blade turbine is designed for a specific 
work output at the mean radius, the average work output will be a maxi- 
mum of approximately 2 percent less than the value based on the mean 
radius . 

Rotor turning angle. - Another factor in turbine design is the 
rotor-blade turning angle. In free-vortex turbines, the turning angle 
is maximum at the hub where the relative rotor -entrance Mach numbers 
and the reaction factor are most critical. The rotor -blade -hub turning 
angle is plotted against hub -tip -radius ratio for the nontwisted-rotor - 
blade and the free-vortex turbines in figure 10. In every case, the 
turning angle, and thus the blade loading, for the nontwisted-rotor - 
blade turbine was less than the turning angle for the free-vortex tur- 
bine, so that it is probable that the use of higher relative rotor-hub - 
entrance Mach numbers and lower reaction factors is possible with 
nontwisted-rotor -blade turbines than with corresponding free-vortex 
turbines. 

Rotor -exit tangential -velocity head. - For nontwisted-rotor -blade 
turbines, it is impossible to obtain zero values of exit tangential 
velocity at all radii. The average specific rotor-exit tangential- 
velocity head (that is, the mass -averaged value of V u ^ 3 2 /2gjj is 

plotted against hub -tip -radius ratio in figure 11 for the cases con- 
sidered. Charging this velocity head as a loss to the turbine is 
equivalent to a maximum decrease of efficiency of the nontwisted-rotor - 
blade turbine of approximately 1 percent compared with the free-vortex 
turbines. 
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Radial Equilibrium Accounting for Streamline Curvature 

In order to obtain an indication of the effects of the radial 
shifts of weight flow on the design parameters , as discussed in the 
preceding section, the analysis accounting for streamline curvature in 
the axial-radial plane was applied to a nontwisted-rotor -blade and a 
free-vortex turbine. The following design conditions chosen were those 
of a typical contemporary turbine: 


Total -pres sure ratio . . . ^ 2.28 

Turbine-inlet temperature, R 2000 

Turbine internal efficiency 0.85 

Weight flow, lb/sec 75.5 

Turbine speed, rpm 11,500 

Tip radius, in 13.0 

Hub radius, in 9.0 

Hub -tip -radius ratio 0.693 

Stator aspect ratio (based on axial chord) 2.0 

Rotor aspect ratio (based on axial chord) 1.78 

Axial clearance between stator and rotor, in 0.72 


The results of this analysis are presented in figures 12 to 15, 
together with the results of a simplified-radial -equilibrium analysis 
based on the same design conditions. The amount of radial shift encoun- 
tered in the axial clearance space midway between the rotor and the 
stator is plotted against radius ratio r/r m in figure 12. It may be 
seen that the deviation of the streamlines is toward the tip for the 
free-vortex turbine and is toward the hub for the nontwisted-rotor -blade 
turbine. In both cases the deviations are larger for the simplified- 
radial-equilibrium analysis. In the free-vortex turbines , however, the 
radial shift of flow is smaller than in the nontwisted-rotor -blade tur- 
bines, especially in the case where the curvature of the streamlines is 
considered. The stator-exit-angle variation for the nontwisted-rotor- 
blade turbine is presented in figure 13 for both the simplified-radial- 
equilibrium analysis and the analysis accounting for streamline curva- 
ture . It may be seen that the variation in stator -exit angle is less 
than 1°. A radial variation of circulation across the stator exists 
with this type of stator. The effect of this variation of circulation 
on turbine performance is unknown. It is noted, however, that the 
radial gradient of circulation was less than the circulation gradient 
of a wheel-type stator where V u = Kr for the same mean- radius condi- 
tions and the wheel-type stator has been used in current compressor 
designs with no decrease in efficiency. 

Relative rotor-entranc e Mach number. - In the discussion of rel- 
ative rotor-entrance Mach number for the simplified-radial-equilibrium 
analysis, it was found that the relative Mach number at the hub was 
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slightly higher for the nontwisted-rotor -blade turbine than for the 
free-vortex turbine. This difference was also found for the typical 
turbine design shown in figure 14 , where the relative rotor -entrance 
Mach number is plotted against radius ratio r/r m for the simplified- 
radial -equilibrium analysis and for the analysis accounting for stream- 
line curvature . On the basis of the analysis accounting for the cur- 
vature of the streamlines, however, the radial pressure variation has 
changed so that the velocities at the hub for the free-vortex turbine 
increased and the velocities at the hub for the nontwisted-rotor -blade 
turbine decreased. The result was an increase in Mach number at the 
hub of 0.05 for the free-vortex turbine and a decrease of 0.08 for the 
nontwisted-rotor -blade turbine. Thus, the comparison of relative 
rotor -entrance Mach number based on the simplified-radial-equilibrium 
analysis may be considered conservative. 

Reaction factor. - The effect of streamline curvature on the 
reaction factor based on a simplified-radial-equilibrium analysis is 
presented in figure 15, where the reaction factor \|r is plotted against 
radius ratio r/r m . On the basis of the analysis accounting for the 
streamline curvature, the reaction factor at the hub of the turbine 
decreased from -0.21 to -0.26 for the free-vortex turbine and increased 
from -0.11 to 0.03 for the nontwisted-rotor -blade turbine. It may be 
concluded that the comparison of reaction factor may also be considered 
conservative. 

The other design parameters of rotor turning angle, specific work 
output, and rotor -exit tangential velocity were also investigated and 
it was found that the variation from the results of the simplified- 
radial-equilibrium analysis due to streamline curvature was negligible. 


CONCLUSION 

A comparison of the turbine aerodynamic characteristics (relative 
rotor -entrance Mach number, reaction factor, rotor turning angle, radial 
variation of specific work output, and specific rotor-exit tangential- 
velocity head) was made on the basis of the simplified-radial- 
equilibrium analysis for nontwisted-rotor -blade and free-vortex tur- 
bines . Further analysis was made of a typical contemporary turbine 
approximating the forces due to the curvature of the streamlines in an 
axial-radial plane to obtain an indication of the effects of the radial 
shifts of weight flow on the results obtained by the simplified-radial- 
equilibrium analysis. 
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It was found that the aerodynamic characteristics of nontwisted- 
rotor -blade turbine are approximately those of a free-vortex turbine 
designed for similar service for values of hub -tip -radius ratios that 
are used in current turbines. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, January 4, 1951. 
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APPENDIX A 
SYMBOLS 

The following symbols are used in this report: 

A annulus area, sq ft 

a cr velocity of sound, based on critical temperature, ft/sec 

b axial distance from point of inflection to point of maximum 
deflection of streamlines, ft 

C radius of curvature of streamline in axial -radial plane, ft 
c p specific heat, ft-lb/(lb) (°R) 

g acceleration due to gravity, ft/sec 2 

h specific enthalpy, ft/lb 

h specific-work ratio 

J mechanical equivalent of heat, ft-lb/Btu 

K function of (g, Cp, and y) 

M Mach number 

p absolute pressure, lb/sq ft 

R gas constant, ft-lb/(°R) (lb) 

r radius , ft 

T absolute temperature, °R 

U blade velocity, ft/sec 

V absolute gas velocity, ft/sec 

W relative gas velocity, ft/sec 

w gas weight flow, lb/sec 

w gas -weight -flow ratio 


16 


NACA El 2365 


x axial distance, ft 

a angle of absolute velocity with, tangential direction, deg 

p angle of relative velocity with tangential direction, deg 

y ratio of specific heats 

A maximum deflection of streamlines in radial direction, ft 
turbine efficiency based on stagnation conditions 

p density, lb/cu ft 

cp angle between streamline and axial direction, deg 

\|r reaction factor 

cd angular velocity, radians/sec 

Subscripts: 

0 reference radius 

1 stator entrance 

2 stator exit, rotor entrance 

3 rotor exit 

h turbine hub 

m mean radius 

r radial 

t turbine tip 

u tangential 

v relative 

x axial 

Superscript: 

1 stagnation state 
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APPENDIX B 

RADIAL VARIATION OF VELOCITY AT ROTOR ENTRANCE FOR 
SIMPLIFIED -RADIAL-EQUILIBRIUM ANALYSIS 

The following analysis of the flow conditions at the rotor entrance 
is made with the specification that the relative rotor -entrance angle is 
constant along the radius. It will he assumed that simplified radial 
equilibrium exists in the axial-radial plane; that is , the forces due 
to the curvature of the streamlines are zero. Axial symmetry, constant 
total temperature along the radius, and constant entropy along the 
radius will also be assumed. 


From vector -diagram considerations (fig. l) , 

v x V x 

tan 13 " TvpuT = (v„-i?)‘ 


(Bl) 


u 


u 


and differentiation of equation (Bl) with respect to the radius gives 


dv x + R ( dV u ) 

_ = tan p 


(B2) 


The energy equation for adiabatic compressible flow may be written as 

2 „ 2 

+ h = h' (B3) 


V “ V 

U X 

2g + 2g 


and differentiation of equation (B3) with respect to the radius, with 
dh' /dr = 0, gives 


V u d ^u ^x dV x + dh _ q 

~g dr g dr dr 


(B4) 


For simplified radial equilibrium^ 

V 2 
v u 

gr 


dp 

p dr 


(B5) 


With constant entropy along the radius , 


dh 

dr 


p dr 


(B6) 
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Equations (Bl) , (B2), and (B4) combine to yield 

2 

1 

g 


^ H + 12sL 1 (v u ^r) & - o>) + 
g dr g u \dr / 


dh 

dr 


= 0 


and equations (B5) and (B 6 ) combine to yield 

v 2 
dh V u 


dr gr 


The substitution of V u 2 /gr for dh/dr in equation (B7) yields 


V u tan p 


h — (V -cur) tan B i - 5 — 

g dr g v u ' p \dr 


dV. 


u I 'u . 

a>) + = 0 

gr 


Simplifying equation (B9) and making the substitution for or — U 
yields 

2 

dV^ V u V u - » tan 2 P + ^ tan 2 P 
r v u - U tan 2 0 + V u tan 2 0 


(B7) 


(B 8 ) 


(B9) 


(BIO) 


When the numerator and the denominator of the right side are divided 
by V u , equation (BIO) may be written as 

t 2 


1 - JL tan 2 p + Hi tan 2 p 
dV u = V u Vu 2 

r tt o p 

1 - — tan^ 0 + tan^ 0 
v u 


dr 


Differentiating V u /U with respect to r yields 

/V. 


dr 


iVi^ 

2 o U dr 


(Bll) 


(B12) 


Substituting equation (Bll) for dVy/dr in equation (B12) and 
dividing the numerator and the denominator by r Q yields: 
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-ft). 


ft 


r 

r 


U 


\3 C 


1 - — tan p + — ^ tan P 


¥ u 

If 


u 


*u 


1 - =^- tan 2 p + tan 2 p 
v u 


(B13) 


The limits -will be assigned as V u /u = 1 to V u /U and r/r Q = 1 to 
r / r o- 

Integrating equation (B13) results in 

v2 


= e 


1/1 + tan^ p\ , 

\ log • 


M _ o( V uV tan 2 p j + tan 2 p 
.U / \U J\2 + tan 2 p / 2 + tan' 1 p 


V 


tan 2 p 


tan 2 p 


V u tan 2 p 




Vff t — 2 1 + tan — P / arc tan 

1 tan 2 B A _ tan 2 P \ 

2 + tan 2 p \ 2 + tan 2 p/ 


2 + tan 2 B 


tan 2 6 /]_ _ tan 2 6 V 

2 + tan 2 p V 2 + tan 2 p/ > 


u 

U 


1 

(B14) 


Substituting the l imi ts in the right member results in 


r 

r^ 


\cos^ p + l y 


log* 


wf + ^ (t - ^ 


= e 


sin 2p , 
— arc tan 

tsfz 


tan p 

tjz 


(* - v 


( 1 ) 
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APPENDIX C 

CONSTRUCTION OF FIGURE 3 (ROTOR ENTRANCE) 

The chart of figure 3 was constructed as follows: An increment 

of weight flow through an increment of annular area may he expressed 
as 

dw = pV x 2jtr dr (Cl) 

The ratio of static to total density may he expressed in terms of the 
critical velocity as 




, N 2l 


n r-i 

v \ 

— 

1 " r+i 

V a cr/_ 


1 


Substituting for p in equation (Cl) yields 



r- 1 i 

f v f 

dw = p 1 

1 - r+i 

V a cr /_ 


1 

Fi 


V x 2jtr dr 


(C2) 


(C3) 


Inasmuch as 


a cr 2 = gC I ,T ' 


V x = (V u -U) tan P 


and 


V = V u 2 + (V u -U) 2 tan 2 p 
equation (C3) may he written as 


dw = p' 


L 


V u 2 + (Vu-U ) 2 tan 2 P 
2gc p T ' 


r-i 


(V u -U) tan p2jtr dr (C4) 
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Substituting a^. for U in equation (C4), dividing by r Q S and making 
the substitution p’ = p'/^'; where R is the universal gas constant, 
yields 


1 



and dividing equation (C5) by R/^/T 7 results in the following expres- 
sion for the weight-flow parameter; 
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For given values of Mo, r/r Q was assigned within the range given 
in figure 2 and values of V u /u were determined for constant values 
of P from figure 2. The equation was then numerically integrated 
within the limits of r/r Q = 1 to r/r Q . The results of the integration 

w a/t** 

are presented in figure 3, where o is plotted as a function of 

p f 7tr 0 ^ 

radius ratio r/r Q for constant values of (3 for x = 1*3. 

/ 

The lines of constant relative -rotor entrance Mach number and 
stator-exit angle are plotted on these figures and were determined as 
follows: 

The relative rotor -entrance Mach number may he expressed as 


v = VrsRT 

From the vector -diagram considerations, 

W 2 = W u 2 + V x 2 = W u 2 + tan 2 3 W u 2 


(C7) 


then 

W 2 = (v u - U) 2 (l + tan 2 p) (C8) 

and from the definition of total temperature, T may be expressed as 


T = T' - ^ (V x 2 +V u 2 ) 


(C9) 


Inasmuch as V x 2 = tan 2 p (V u -U) 2 , equation (C9) may be written as 


T = T* - 


r-i 

2rgR 


[jv u -U) 2 tan 2 p + V u 2 ] 


(CIO) 
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When equations (C8) and (CIO) are substituted in equation (C7) squared, 
the following expression is obtained for My 2 : 




2 


(V u -U) 2 (1 + tan 2 p) 


(Cll) 


rgRT' - [tan 2 p (V u -U) 2 + Vu 2 ] 


Dividing the numerator and the denominator of the right member by U 2 


and making the substitution of U 2 = 

o 

the definition of M 0 yield 



o 


2 and the substitution of 


r+ 1 'ro\ 2 _ 
2 Vr / M_ 



(C12) 


From values of Vu/U determined from figure 2 for assigned values of 
r/r G and P, M w was calculated from equation (C12) . When M w was 
plotted against r/r Q for constant values of p and Mq , lines of 
constant relative rotor-entrance Mach number were obtained, as presented 
in figure 3. 

The variation of stator-exit angle with radius ratio r/r Q was 
determined as follows: 

From vector -diagram considerations, 



or 


1 


(C13) 


u 


tan a 


1 


tan p 
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When constant values of a and. p were assigned within the range 
used in figure 3, Vu/U was calculated using equation (C13). The cor- 
responding value of r/r G was determined from figure 2 and lines of 
constant stator-exit angle were plotted against r/r Q in figure 3. 
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APPENDIX D 


RADIAL VARIATION OF TANGENTIAL VELOCITY AT ROTOR EXIT FOR 


SIMPLIFIED-RADIAL-EQUILIBRIUM ANALYSIS 


For the analysis at the rotor exit, as in the analysis at the 
rotor inlet, it is assumed that simplified radial equilibrium exists 
in the axial -radial plane; that is, equation (B5) applies. 

The energy equation for adiabatic flow along a streamline between 
the rotor entrance and exit may be written as 


Inserting the expression given for work into equation (Dl) and then 
differentiating with respect to r yield the following equation: 



= h ' g - work 


(Dl) 


Work may be expressed as 


Vu,2 U 2 Vu,5 u 5 


(D2) 


g 


g 



dr g dr g dr 



(D3) 


drgdrgdr gdr gdr 


Constant entropy along the radius is again assumed; therefore 
appendix B, 


V u , 3 2 


dr gr 


as in 
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Inasmuch as 


and 


v x,3 = tan ^3 ( v u,3" u 3) 


Wx,3 A „ ( d V u ,3 

-±r = tan Ps vrar 


®3 

dr, 


equation (D3) may be written as 


v u,3 2 v u,3 dv u,3 tan2 h („ _ dV u,3 

+ ~r — aF" + g v u > 3 


II dV u,3 
U 3 dr“ 


dU 3 dU 3 

^U,3 (Jp ^3 dr, 


^'2 u 2 dv u,2 u 3 dv u,3 V u,2 ^2 | V u,3 dU 3 

“dr g" dr“ + "g~ dr “ g dr + g dr 


(D4) 


dh' 


Because 


dr 


dV n 


= 0, solving for 


u,5 


dr 


in equation (D4) yields 


dV, 


u,5 


dr 


v -r au, „ au, dv„ ? au 2 „ du 3 

+ V u,5 d# + tan p 3 V u,3 - U 2 ~ V u,2 dr~ ~ tan P 3 U 5 W 

V u,3 + V u,3 tan " P 3 ' U 3 tan2 P 3 * U 3 


(D5) 


When terms are combined, equation (D5) may be written as 


dV. 


u,5 


v 2 r „ n du 3 

- — + |_ v u,3 + tan 2 P 5 (V U;3 -U 3 )J - U 


dV, 


dUc 


uzl _ v „ ri 


2 dr " u )2 dr 


dr 


(Vu,3" u 3) (1 + tan P3) 


(D6) 
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When it is assumed that Ug = U3, dividing the numerator and the 
denominator of the right member by co results in 


avu,5 

dr 


- Aklf + v u , 3 (1 + tan2 p 3 ) - r 


dVu,2 


dr 


- V, 


u, 2 


U tan^ 


P3 


( 3 ) 


V, 


u,3 

~ 


■lJ\l + tan 2 p 3 


) 
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APPENDIX E 

CONSTRUCTION AND USE OF FIGURE 4 FOR DETERMINATION OF 

ROTOR-EXIT ANGLE 

The chart of figure 4 for determining the rotor-exit angle was 
constructed as follows: From the energy equation 


1 


p 

-1 T-l 

( V \ 2 

r-i 

p ! 

1 -riT 

\ a cr / 



where 


a cr 



and from the exit -velocity vector diagram 


.2 / \t \2 / -it \ 2 

V ) / v x ) | *u 


cr 


cr 


cr 


the following expression is derived: 


pV : 


x 


p’ a cr 


1 - 


--1 /V x 2 + Vu 2 


r+1 


x cr u 


r - 1 


v. 


X 


cr 


(El) 


For given values of ranging from 0 to 1.2 and values of 

a cr 

Vx ranging from 0 to 1.0, the absolute velocity ratio and cor- 


ner 


pVx 


rected weight-flow parameter were calculated for y - 1*30. The 

P a cr 
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results are presented in figure 4, 

with lines of constant values of 

of constant absolute flow angles, 
ing expression: 


V pV x 

where is plotted against — ; 

a C r P a cr 


Vx 

a cr 


and 


Vu 

a cr 


Superimposed are lines 


which were determined from the follow - 


tan“^- a 


Vx/ a cr 
Vu/ a cr 


The procedure for applying figure 4 to determine the rotor-exit 
angle at the mean radius is as follows: With V u ,3, assigned at the 

mean radius, the work was calculated from equation (D2). The stagnation 
temperature and pressure at the rotor exit were then determined from 


and 


T . = T * - 

3 1 c p J 



r 

work ^ 
Cp J rjipT ' p j 


Values of p ' 3 and a cr ^ were determined from the values of T'j 
and p', thus obtained. When it is assumed that the specific weight 
flow at the mean radius is equal to the average specific weight flow, 
the corrected weight-flow parameter may be determined from the following 
expression: 


f— = s£ 

\ P ' a cr/ 3 % p ' 3 a cr f 3 


Then 

value 


(-) 

N a cr /3 

of 33 


is obtained from figure 4 and, when U 3 
may he calculated. 


is known, the 
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APPENDIX F 

RADIAL VARIATION OF VELOCITY AT ROTOR ENTRANCE FOR RADIAL 

EQUILIBRIUM ACCOUNTING FOR STREAMLINE CURVATURE 

The radial-equilibrium equation at a station midway between the 
stator exit and the rotor entrance , when constant entropy along the 
radius is assumed, may be written as (from reference 3, equation (14)) 

dh' _ ^u d (y) v * * * * x dv x v x dv r 
dr “gr dr + g dr ~ g dx 


when is assumed equal to zero, equation (Fl) may be written as 

dr 


V <3.V r dV u 

— V v — - + v„ -s — + v 

r x dx u dr 


dV 

x dF 


x 


= 0 


(F2) 


The radial component of velocity V r = V x tan cp. (See fig. 1.) Dif- 
ferentiating V r with respect to x and substituting the result in 
equation (F2) results in 



2 d( tan cp) 
" v x dx 


dV. 


dV. 


dV. 


- V x tan cp ^ + v u dr + V * dr 


= 0 


(F3) 


It will be assumed that the maximum deflection of the streamlines 

occurs at this station. Therefore, cp equals zero and 

d(tan cp) _ dcp _ 1 equation (F3) may then be written as 
dx dx C 



r 


V x 2 dV u dV x 

FT + V u dF" + v x dF" 


0 


(F4) 


Nontwist ed-rotor -blade turbines. - From vector -diagram considera- 
tions for nontwisted rotor blades, 

v x 

tan p = 77 n = constant 

v u - u 


or 


V u = V x cot p + U 


(F5) 
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Differentiation with respect to the radius gives 


dV u dV x 

aF~ = cot p dF" + 03 


(F6) 


Substituting equations (F5) and (F6) in equation (F4) and sinrpli- 

dV x 

fying result in the following expression for 


dV, 

dr 


(V x cot p + U) 


V 


+ 0 ) (V x cot 3 + U) 


V x + v x cot^ p + U cot p 


(F7) 


Using the identity 


dV 


x 


dr 


SU 


dr 


\ dU 
U dr 


multiplying equation (F7) by r m /u and the numerator and the denominator 
of the right member by l/u, and simplifying yield 


2+4: 


(aT " 

\r m / r m 


1_ 

r m 


cot p + 




m. 


+ (1 + 2 cot p) 


V x P 

— CSC^ P + cot p 


^ (4) 


Free -vortex turbines. - For free -vortex flow. 


V u r = constant 


u 


Differentiating with respect to r and multiplying by — yield 




(F9) 


Substituting equation (F9) in equation (F4) and simplifying result in 


^x _ Vx 
dr C 


( 5 ) 
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Figure 1. - Velocity vector diagram and section of turbine 
through axial-radial plane. 
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- Variation of tangential-velocity ratio with radius ratio for constant rotor- 
entrance angles. (A 17- by 22-in. print of this fig. is attached.) 
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Figure 3. - Variation of weight-flow parameter, stator-exit angle, and relative 
rotor-entrance Mach number with radius ratio for constant relative rotor- 
entrance angle (7 = 1.3). 
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Figure 3. - Continued. Variation of weight-flow parameter, stator-exit angle, 
and relative rotor-entrance Mach number with radius ratio for constant relative 
rotor-entrance angle (y = 1.3). 
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(c) Mq = 0.8. 

Figure 3. - Continued. Variation of weight-flow parameter, stator-exit angle, 
and relative rotor-entrance Mach number with radius ratio for constant relative 
rotor-entrance angle (y = 1.3). 
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(d) M 0 = 0.7. 

Figure 3. - Continued. Variation of weight-flow parameter, stator-exit angle, 
and relative rotor-entrance Mach number with radius ratio for constant relative 
rotor-entrance angle ( y = 1.3). 
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Figure 3. - Continued. Variation of weight-flow parameter, stator-exit angle, 
and relative rotor-entrance Mach number with radius ratio for constant relative 
rotor-entrance angle (7 = 1 . 3 ). 
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Figure 3. - Concluded. Variation of weight-flow parameter, stator-exit angle, 
and relative rotor-entrance Mach number with radius ratio for constant relative 
rotor-entrance angle (7 = 1.3). 
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Figure 4. - Variation of absolute velocity ratio with corrected weight-flow ratio for 
constant absolute tangential -velocity ratio and absolute axial -velocity ratio (*r = 1.3). 
(A 17- by 22-in. print of this fig. is attached. ) 
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Figure 5. - Variation of weight-flow ratio 
with hub-tip-radius ratio for nontwisted- 
rotor-blade and free-vortex turbines for 
same mean-radius conditions. 
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Figure 6. - Variation of hub relative 

rotor-entrance Mach number with hub-tip- 
radius ratio for nontwisted-rotor-blade 
and free-vortex turbines for equal 
weight flow per unit annulus area. 


Reaction factor at hub, 


44 


NACA TN 2365 



Figure 7. - Variation of hub reaction 
factor with hub-tip-radius ratio for 
nontwisted-rotor-blade and free-vortex 
turbines for equal weight flow per unit 
annulus area . 
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Figure 8. - Variation of ratio of specific work to specific 
work at mean radius with radius ratio for nontwisted rotor- 
blade and free-vortex turbines. 
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Figure 9. - Variation of average-specific- 
work ratio with hub-tip-radius ratio 
for nontwisted-rotor-blade and free- 
vortex turbines for equal weight flow 
per unit annulus area. 
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Figure 10. - Variation of rotor-hub 

turning angle with hub- tip-radius ratio 
for nontwisted-rotor-blade and free- 
vortex turbines for equal weight flow 
per unit annulus area. 
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Figure 11. - Variation of specific rotor- 
exit tangential-velqcity head with hub- 
tip-radius ratio for nontwlsted-rotor- 
blade and free-vortex turbines for 
equal weight flow per unit annulus area. 
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Radius ratio, r/r m 

Figure 12. - Variation of radial deflection of streamlines midway 
between rotor and stator with radius ratio for typical turbine 
example . 
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Figure 13. - Variation of stator-exit angle with radius ratio for 
typical example of nontwisted-rotor-blade turbine. 
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Figure 14. - Variation of relative rotor-entrance Mach number with 
radius ratio for typical turbine examples. 
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